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Abstract. 

We discuss collision of three identical particles and derive scattering selection 
rules from initial to final states of the particles. We use either hyperspherical or 
Jacobian coordinates depending on which one is best suited to describe three different 
configurations of the particles: (1) three free particles, (2) a quasi-bound trimer, or (3) 
C^ I a dimer and a free particle. We summarize quantum numbers conserved during the 

lyj ' collision as well as quantum numbers that are appropriate for a given configuration but 

may change during the scattering process. The total symmetry of the system depends 
£^ I on these quantum numbers. Based on the selection rules, we construct correlation 

(— I ' diagrams between different configurations before and after a collision. In particular, 

^ I we describe a possible fragmentation of the system into one free particle and a dimer, 

which can be used, for example, to identify possible decay products of quasi-stationary 
^SJ ■ three-body states or three-body recombination. 
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1. Introduction 

Quantum systems of three identical particles have attracted a considerable interest 
during the last years, principally due to very recent experimental and theoretical 
advances in the field of ultra-cold degenerate gases. One such example of the recent 
experimental advances is the observation of three-body Efimov states [1] in ultra cold 
gas of cesium [2]. Another motivation to study quantum few-body problems is due to 
recent experiments with interacting atoms and molecules trapped in optical lattices. 
The interaction between a few atoms (or molecules) in optical traps can be changed by 
an appropriately chosen external field. Understanding and controlling the interaction 
between atoms (and molecules) in optical lattice contributes into the progress in the 
domain of quantum information and quantum computing. 

One aspect of the interaction in the three-body system is the selection rules 
determining allowed and forbidden final states of the system after the scattering process. 
For three identical particles, there are three possible configurations before, during and 
after a collision. Namely, (1) the three particles are far from each other and do not 
interact at all, (2) two of the three particles are close and interact with each other, but 
the third particle is situated far from them and therefore can be considered as free, 
and (3) the three particles are close to each other such that the interaction between 
them cannot be neglected. For each of these configurations, one can assign a set of 
appropriate quantum numbers. Before and after a collision, the configuration may 
change. Therefore, the set of quantum numbers may also be different before and after 
the collision. However, there are relationships between the sets corresponding to different 
configurations. In addition, some quantum numbers are the same in the three sets and 
are therefore conserved during the collision. Such correlations between the three sets 
of quantum numbers are the consequence of the conservation of symmetry of the total 
wave function of the system. The symmetry of rotational and vibrational wave functions 
of the system of three indistinguishable particles has been studied previously by several 
authors [21 SI [5l El [7l E] • However, there was no systematic description of the correlation 
between quantum numbers appropriate for all the three possible configurations. In this 
article, we discuss the three sets of appropriate quantum numbers and derive the rules 
describing the correlations between different quantum numbers in the three sets. Since 
there exist some selection rules regarding the change in the quantum numbers, we derive 
and discuss them in detail. 

In section [21 we will give a brief overview of symmetry properties of the three 
identical particles and general rules for construction of wave functions. Then, in section 
131 we present wave functions and the corresponding quantum numbers appropriate for 
each of the three configurations of the system. In section HI we derive the scattering 
selection rules and represent them by diagrams of correlation between the three sets 
of quantum numbers, similar to the case of diatomic molecules [9j. Finally, section [5] 
contains our conclusions. 
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2. System and symmetry 

The group S3 of permutations of three identical particles is isomorphic to the point group 
of symmetry C^^. Geometrically, the system with the C^^ symmetry can be represented 
by a flat equilateral triangle containing the six following elements: three reflexions cr„ 
through three vertical planes perpendicular to the plane of the triangle, two rotations C3 
by 27r/3 in the plane of the triangle and the identity operator E [TOj [11]. Adding to S3 
the operator of inversion /, the new group 6*3®/ becomes isomorphic to the point group 
D^h, which itself can be represented as D^h = Csv^cTh, where (7^ is the reflection through 
the symmetry plane. Relative motion of three structureless particles is characterized by 
the C^y group. If rotation of the whole system in a space-fixed coordinate frame needs 
to be accounted for, then, the symmetry of the system is represented by the D^h group. 
In our discussion we will be using the same notations for elements g of the group D^h 
as in [11]; the 12 elements g are E, (12), (13), (23), (123), (132) without inversion plus 
the same elements with the inversion E*. 

Each eigenstate of the Hamiltonian of the system is transformed in the D^h group 
according to one of the six D^h irreducible representations: A[,A'l,A2,A2,E', or E". 
The representations E' and E" are two-dimensional. For convenience, we define bases 
in the E' and E" representations. Namely, in the two-dimensional space of E' we use 
two orthogonal basis states, E'j^ and E'_, whose properties are such that 

(123)E± = e'^^E'^ (1) 

(12)^1 = E^, (2) 

where u = 27r/3. The basis states turn into each other through operator (12). The 
basis states, E^ and E'_!_ in the E" representation are defined in the same way. 

Having defined the bases in E' and E", we can represent an arbitrary state 
i/j of the system as a superposition of eight functions, each transforming as 
A[,A'l,A2,A2,E'_^,E'_,E1, or E'^. Each term in the superposition can formally be 
obtained by projecting the state i/j on the corresponding function. The projector 
operator on the basis state T is represented by the following sum 

Pr(^)= Yl ^Jg^- (3) 

Here, the sum is taken over all 12 elements g of the D^h group. This formula is just the 
regular projector on a one-dimensional irreducible representation F of D^^. Therefore, 
in the above equation, x^ are the well-known characters [TI] if F is an one- dimensional 
irreducible representation. Coefficients Xg for the states E'^ or E'^ are given in table [H 

For simplicity, we will only consider the coordinate part of total wave function of 
the system. The electronic and nuclear spin statistics of three identical particles can be 
described using a similar approach, which is briefly discussed in our conclusion. 

We are interested in selection rules in the collision process, which starts with a given 
state xp of the system. The state ip is usually specified by a set of quantum numbers 
that are appropriate for the initial configuration of the system. These quantum numbers 
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Table 1. Coefficients Xg in the projector formula ([3]) for _Ej_ and E'^ basis states. 
Coefficients x^ for one-dimensional irreducible representions of the D^h group are 
equal to the corresponding characters [11]. Operators E and E* are the identity and 
inversion operators correspondingly. 



are in general not conserved during the scattering process. However, they are sufficient 
to decompose the state into the sum of irreducible representation of the group D^h- A 
given irreducible representation is conserved during the collision and can be viewed as a 
good quantum number, allowing the characterization of the possible final states of the 
system. Therefore, decomposition of initial states into irreducible representations of D^h 
allows us to derive the selection rules, i.e. to point out which final states (characterized 
by other set of quantum numbers) are not accessible from that given initial state; if 
the operator Pp acting on ip gives identically 0, it means that final states transforming 
according to F are not accessible from the initial state ip. 

3. Configurations and adapted wave functions 

Before deriving the selection rules, we derive relations between appropriate quantum 
numbers at different configurations of three identical particles and irreducible 
representations that the quantum numbers correspond to. These quantum numbers 
describe entirely symmetry properties of the spatial wavefunction. 



3.1. Short distances 

We consider that the three particles are situated at short distances from each other if the 
potential of interaction in each pair of the particles is much larger than the rotational 
energy of the relative motion. We define the short range region in this way to be able 
to separate vibrational and rotational motion of the system and use vibrational and 
rotational quantum numbers of the three particles. Therefore, the total wave function 
is written as a product of rotational TZ and vibrational (p factors: 

V' = 7^(a,/3,7)0(p,^,<^), (4) 

where a,/3,7 are the Euler angles defining orientation of the plane of the three 
particles in the space-fixed coordinate frame; p, 6, and ip represent three independent 
vibrational coordinates defining relative positions of the three particles. Here, we use 
the hyperspherical coordinates [I2] to describe the relative motion. 



Correlation diagrams in collisions of three identical particles 



The rotational wave function of the D^^ system is described by the symmetric top 
eigenstates [TO] : 
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where J is the total angular momentum of the system, M and K are respectively its 
projections on the space-fixed and molecular axes. These numbers are constants of 
motion of the symmetric top and sufficient to define completely the rotational state 
of the system. The correspondence between the quantum numbers J, M, K and the 
irreducible representations of such rotational functions is summarized in table [21 which 
is obtained from the symmetry properties of the Wigner functions D^^ ^ (a, /3, 7) [TT|, [ 
in the D^f^ group and from the coefficients Xg of ©• 
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Table 2. Allowed quantum numbers for the rotational wave functfon 7?.(a,/9, 7). 

In the rigid rotor approximation, the vibrational wave function depending on three 
coordinates is transformed according to the irreducible representations of the C^^ group. 
If such vibrational states are obtained numerically, their irreducible representations 
are determined by their symmetry properties. If the normal mode approximation is 
used for the vibrational states, one can assign to each vibrational state the label of 
an irreducible representation based on three vibrational quantum numbers fi,f2 and 
I2 describing the three-dimensional harmonic oscillator of the C^^ symmetry. We use 
notation |fi,f2^) to specify the corresponding wave function. The quantum number 
Vi describes the symmetric stretch mode, which preserves the equilateral configuration 
of the system. The value of V2 indicates the number of asymmetric radial vibrational 
quanta. The quantum number I2 indicates the number of quanta in the vibrational 
angular momentum mode and can only take values — f2, —V2 + 2 . . .^2 — 2,f2. The 
transformation properties of these vibrational functions are well known (see, for example, 
[m [151 [IE] ) and depend only on the quantum number I2 (see table [3]). 

An example of the system where normal mode approximation works well for several 
low-lying vibrational states is the Hj]" ion. For instance, the Hjj' wave functions and the 
corresponding normal mode quantum numbers are shown and discussed in [16]. The 
bound state component of three-body resonances could also possibly be represented 
using the normal mode approximation. However, the approximation cannot be used 
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Table 3. The table summarizes the C^^ symmetry properies of the vibrational wave 
functions obtained in the normal mode approximation. 



for the continuum component of the resonances and for scattering states of the three 
particles. They have to be calculated numerically. 



3.2. Large distances 

We consider that the particles are at large distances if at least one particle is far from the 
two others such that the interaction of the particle with the two others can be neglected. 
There are two distinct configurations at large distances: three free particles or one free 
particle and two interacting particles. For simplicity, we call the interacting pair dimer. 
Correspondingly, the Hamiltonian of the system at large distances can be separated in 
two or three non-coupled Hamiltonians describing the free particle(s) and the dimer. 

3.2.1. Three free particles The total Hamiltonian depends on six coordinates (the 
motion of the center of mass is separated): three hyperspherical coordinates and the 
three Euler angles [5], |12] : 
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p^ dp \^ dp J p 

The grand angular momentum operator A depends on the three Euler angles and two 
hyperangles 9 and yj. The symmetry of the eigenfunctions of the Hamiltonian is entirely 
determined by eigenfunctions F^''"''^'^^'^{a, P, 7, (p, 6) of A^ 

The functions _p^'"'-^'*'^'/^ are characterized by five quantum numbers A, t;, J,M,p, each 
corresponding to an operator [5j: 

(i) the grand angular momentum squared A^ with eigenvalues A(A + 4) where A = 
0,1,2,..; 

(ii) the operator S = i2d/dip, which determines the property of F^'"'"^'^^'^ with 
respect to binary and cyclic permutations of the particles and has eigenvalues 

t; = A,A-2,...,-A + 2, -A; 

(iii) the square of angular momentum J^ with eigenvalues J(J + 1) where J < A; 

(iv) the projection Jz of angular momentum on a space fixed axis with eigenvalues 

M = -J, -J + 1,..., J-1, J; 



Correlation diagrams in collisions of three identical particles 7 

(v) the operator lifting the degeneracy for hnearly independent solutions still existing 
for given values of A, v, J, M with the quantum number fi. 

Functions F are written as the finite sum [5]: 



nA,u,J,M,/i/ 



a, P, 7,^,0) 
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Since the hyperangle 9 is unchanged by permutations, the explicit description of the 
function g^^ ^{0) is not important for our discussion. Using the symmetry of the Wigner 
functions and the transformation of the hyperangle (/9, one obtains transformations of 
the function F under the operations of D^h. 

-p* pA,v,J,M,fi __ f T\ApA,v, J, M,^i /q\ 

(12)F^'^''^'^^'^ = (—lyF^~^'^^'^'^^ (10) 

f-i no") rpA,v,J,M,fj, —iujv TpA,v,J,M,fi (^^\ 

The symmetry properties of functions F^'""'"^'^^'^ are summarized in table HI 
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Table 4. Allowed quantum numbers for three free particles 

Since the three particles do not interact with each other, instead of hyperspherical 
coordinates, one can also use spherical coordinates for each particle separately fl, r^, r^. 
Then, the total wave function is written as product 

V'(n, rl, r~^3) = 0a(n)(/)6(r^)0c(rl) , (12) 

where indices a, b, and c specify one-particle states. To obtain definite irreducible 
representations of D^h from the product above one has to apply (IHl). The allowed 
irreducible representations depend on the quantum numbers a, b, and c: If angular 
momenta la, h, h are defined for each state a, b, and c, the parity of the whole system 
is definite and given by (^—ly^+^b+ic _ jf g^^j ^^q q^ ^j^g three states a,b, and c are not 
equal to each other, all three irreducible representations, Ai, A2, or E of the allowed 
parity are possible (the sum of ([3]) does not vanish). If two of the three indices a, b, and 
c are equal to each other, only Ai and E are possible. Finally, if all three indices are 
the same, only Ai irreducible representation is allowed. 
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3.2.2. One dimer and one free particle We start with the non-symmetrized case, let 
say, particles 1 and 2 form a dimer, particle 3 is far from the dimer. The interaction 
between the dimer and particle 3 is negligible. Here, we use the mass-scaled Jacobi 
coordinates defined as r = (3/4)-^/^(r2 — "H,) and R = (4/3)^/^(^^i^ — f^) [3]. In these 
coordinates, the kinetic energy operator has the simple form: 

T = -^(A,^+A^), (13) 

with the 3-body reduced mass /x = m/-\/3. Thus, the Hamiltonian is 

H = -^A^ + Hi,, (14) 

where if 12 describes internal motion of the dimer. The conserved quantities for this case 
are the two angular momenta jr and Jr and their projections m,, and Mn associated 
with r and R correspondingly, the vibrational state v of the dimer and the kinetic energy 
e associated with motion along R. The wave function of the system (excluding the factor 
describing motion of center of mass) has the form 

where Y^-^m(^) is a spherical harmonic depending on two angles collectively called Q, 
(pvjrir) is the vibrational state of the dimer, uj^^e{,R) is the i?-radial part corresponding 
to the free motion, i.e. the usual Bessel function. For us it is more convenient to 
construct the state with a definite total angular momentum J = jr + Jr from the states 
of (US]) as in [3] 

where the bipolar spherical harmonic X-^' j (fi^, ^r) is 

rrir ,Mji 

The wave function of flT6|) does not belong to a definite irreducible representation. 
Again, ([3]) is used to project the wave function on any of the one-dimensional irreducible 
reppresentations or to basis states E'^ and E'^ of D^h- For the present case, it is more 
convenient to rewrite ([3]) in the form where operators (123) and (132) are written as 
products (23)(12) and (31)(12): 

Pr = E + x[2(12) + xL(23) + xLx[2(23)(12) + xL(31) + xLx[2(31)(12) 
+p.t.* = (E + xL(23) + X3i(31)) (E + x[2(12)) (E + Xe^E*) , (18) 

where '+p.t.*' means that we add all previous terms multiplied with the inversion 
operator E*. The effect of operator E* and (12) is easily evaluated: The only factor in 
flT6|) that changes under E* and (12) is the bipolar harmonic. The operator E* acting 
on the VLr part in ( fTTl) gives the same factor multiplied with (— l)"'^, acting on the Vtr 
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part gives an additional factor (— I)-'''. The operator (12) changes only the Vlr factor, 
giving the factor (—I)-''': 

E*<;i = (-1)''--^^'^<;:J, (19) 

Therefore, we can rewrite (ITSl) as: 

PI JM,e 

(1 + xL(23) + xL(31)) (1 + i-iy^x\2) (1 + (-1)'^+^'X^0 <;,.« (20) 

The above combination is identically zero if either the second or the third factor is zero. 
Thus, we obtain the combinations of Jr and jr allowed for states transforming as F. 
The results are summarized in table [51 
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Table 5. Allowed quantum numbers for the dimcr+free-particle configuration. The 
trivial A'l state with jr = J^, = is not included this table. 



4. Correlations between configurations 

The correlations between quantum numbers appropriate at different configurations of 
the three particles are derived from conservation of the symmetry, i.e. irreducible 
representation of the state of three particles. The possible correlations for several values 
of quantum numbers are summarized in figures [1] and [2] for total angular momentum 
J = and J = 1. Since the symmetry depends only on the quantum numbers K and 
I2 at small distances, on A and v for the non-interacting particles, and on jV and Jr for 
the dimer+particle configuration, only these quantum numbers are specified in figures 
[Hand El 

With help of figures [31 and [H we give an example of how the selection rules can 
be applied. Figure [3l shows several adiabatic hyperspherical curves for a three-body 
model system. Such curves are obtained diagonalizing the rovibrational Hamiltonian 
in the space of the hyperangles, treating the hyper-radius as a continuous parameter. 
Rotational degrees of freedom (Euler angles) are also quantized. The particular example 
we consider is the actual Hes potential |T7] multiplied uniformly with an factor to 
obtain a deeper potential well for the trimer. The origin of the energy corresponds to 
the energy of the three particles at infinite separation from each other. For simplicity 
of the discussion, we consider that the total angular momentum J of the system is 
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Figure 1. Correlation diagram between the dinier+free-particle and three-free- 
particles configurations. The quantum numbers used in this diagram are defined in 
section [31 



zero. Therefore, the symmetry of the rotational part of the wave function is A[. 
In figure [3l asymptotically-large values of the hyper-radius correspond to three-free- 
particle configuration if the energy of the system is positive and to the dimer-|-particle 
configuration if the total energy is negative. The lowest possible energy in the asymptotic 
region is the energy of the dimer in the ground rovibrational state. For the considered 
model system, the dimer has only one possible vibrational state v = and four rotational 
levels jr = 0, 1, 2, and 3. 

The selection rules could be applied to collisions of three particles or to the decay of 
a quasi-bound three-body state (half-collision). As an example, we consider the decay 
of a quasi-bound state with the total energy larger than zero. After the decay, the 
system may end up in the three-free-particles configuration or in the dimer-|-particle 
configuration. For example, if the initial state has symmetry A[, the dimer -|-particle 
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Figure 2. Correlation diagram between the dimer+free-particle and three-free- 
particles configurations. The quantum numbers used in this diagram are defined in 
section [31 



configuration is possible only if jr is an even number. This is clearly seen in figure [3l 
adiabatic curves of the A[ symmetry correlate at large hyper-radii only with dissociation 
thresholds with jr = and 2, which is consistant with the correlation diagram shown in 
figure [2j Similar analysis can be done for the E' initial state. However, E' short-range 
states can end up with all possible values of jr- There is no quasi-bound states of the 
A2 symmetry for the considered three-body potential. However, if we assume having 
the three free particles collide in an A2 initial state, after such a collision only dimers 
with odd angular momenta jr can be produced: A2 curves in figure [3] correlate only 
with odd jr at large distances. Such initial arrangement of the three colliding identical 
particles can be produced experimentally, for example, in ultra-cold gas of fully polarized 
fermionic atoms, such as ^°K. 

As mentioned above, the irreducible representation of wave functions in 
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Figure 3. (Color online) Adiabatic curves for the model three-body potential 
discussed in the text. The total angular momentum J = 0. All coordinates except 
the hyper-radius are quantized. Since the irreducible representational is determined 
by the symmetry of the wave functions with respect to the (hyper-) angles ip and 7, 
each curve is labeled with the index of an irreducible representation. For J = 0, only 
positive parity {A[,A2, or E') is allowed. 



hyperspherical coordinates is determined by the behavior of the functions along the 
hyper-angle if and Euler angle 7. In the considered example, the dependence on 7 is 
trivial, since J = 0. Figure H] demonstrates the dependence of the wave functions on 6 
and ip obtained numerically for three qualitatively different hyper-radii: (i) short hyper- 
radius, where the system behaves approximately as a rigid rotor, (ii) the intermediate 
region, and (iii) large hyper-radius, where one can approximate the system by three non- 
interacting particles or the dimer-|-particle configuration. In the intermediate region no 
approximation may be applied. As it is clear from the figure, the symmetry of the wave 
functions is determined by a single variable (p. The binary permutations (12), (23) and 
(13) are equivalent, respectively, to reflexions of the functions about the three axes going 
through the center of each plot in figure H] at angles ip = 7r/6, ip = 57r/6 and p} = — vr/2, 
whereas the cyclic permutations (123) and (132) are respectively equivalent to change 
of the hyperangle from p> to ip + vr/S and (p — vr/S. 



5. Summary and conclusions 

We have considered the sets of quantum numbers of three indistinguishable particles 
for three different possible configurations: three non-interacting particles, the pair of 
two interacting particles and the third one being far from the pair, and finally, three 
interacting particles. Some of the quantum numbers from the three sets are the same in 
the three configurations, others are appropriate only for the given configuration. During 
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Figure 4. (Color online) Hyperangular wave functions of three irreducible 
representations Ai,A2, and E. E is represented by its two components Ea = 
Real(i?+), and E^ — Im(£'+). Here, the parity is defined (positive) but we do not 
specify it here because it is controlled by 7, that cannot be shown. In the left column 
the hyper-radius is « 9 a.u., in the middle column it is « 20 a.u., and in the right 
column it is « 30 a.u. The wave functions arc shown as contour plots depending on 
9 and ip: The three upper functions arc shown for only two values of contour cuts: 
0.5 and 0.75; the other functions as shown for values ±0.5 and ±1. The coordinates 9 
and ip in this graph are arranged in the polar-coordinate-like way: 9 plays the role of 
radius and 4> plays the role of cyclic polar angle. 



the process of collision of the three indistinguishable particles, all three configurations 
are involved. The quantum numbers that are the same for the three configurations, are 
conserved during the collision, others are transformed. We analyzed the correspondence 
between quantum numbers from the three sets and derived certain selection rules 
between them. The selection rules for several values of the quantum numbers are 
schematically presented in figures [H and [2] as correlation diagrams similar to the 
correlation diagrams derived for diatomic molecules [9] . If needed, correlation diagrams 
for other values of the quantum numbers can be easily obtained from tables E], [3l HJ 
and [5l The derived correlation diagrams can be useful in analysis of experiments with 
ultra-cold degenerate atomic gases studying such processes as three-body recombination, 
formation and decay of three-body resonances [21 [181 [13 EO, EI], or thermalization of 
the degenerate gas. 

Concluding, we would like to discuss briefly how the spin statistics can be included 
in the selection rules. Treating the nuclear and electronic spin parts of the total wave 
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function separately from the rovibrational part, we neglected the interaction of the 
rovibrational motion with (electronic and/or nuclear) spins of the particles. Hence, 
we assumed that the spin Hg Hamiltonian is separable from the the rovibrational 
Hamiltonian and the total wave function of the system is constructed as a product 
of the spin and spatial factors. The symmetry group of Hg is the subgroup C^v of D^h. 
Correspondingly, spin eigenstates transform in D^h according to one of the A[, A2, or E' 
irreducible representations. The effect of the total spin of the system on the symmetry 
of total wave functions can be accounted using well-known tables of products of the D^h 
irreducible representations. An example of relationships between nuclear spin quantum 
numbers and irreducible representations for the system of three identical particles is 
given in ^Q\ for the case of Hg" (spin-^ particles) and D^ (spin-1 particles). 
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